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We examine the behavior of the Kolmogorov constants C2, and C k i, which are, 
respectively, the prefactors of the second order longitudinal structure function, the three 
dimensional and one-dimensional longitudinal energy spectrum in the inertial range. We 
show that their ratios, CijC k \ and Ck/Cki, exhibit clear dependence on the micro-scale 
Reynolds number R\, implying that they cannot all be independent of R\. In particular, 
it is found that {C kl /C 2 - 0.25) = 1.95-R^ 68 . The study further reveals that the widely- 
used relation C2 = 4.02Cfci holds only asymptotically when R\ > 10 5 . It is also found 
that C2 has much stronger independence than either C k , or C k i if the latter indeed has 
a systematic dependence on R\. We further show that the variable dependence on R\ 
of these three numbers can be attributed to the difference of the inertial range in real- 
and wavenumber-space, with inertial range in real-space known to be much shorter than 
that in wavenumber space. 
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1. Introduction 



The idea of homogeneous and isotropic turbulence (HIT) (|Kolmogorovlll94ll ) enables 
one to focus on the essential physics of small-scale turbulent properties in the simplest 
possible case and use it as a first step to understand more c omplicated turbulence prob- 



lems The Kolmogorov 1941 model (K41) (jKolmogorovl 1941 ) for HIT provides predictions 



of turbulence properties that agree well with experiments when the effect of intcrmittency 
is negligible. One example is the "two-thirds law" for second-order velocity structure func- 
tions and its counterpart in wavenumber space, the "five-thirds law" for energy spectra, 
all expected to hold in the so-called inertial range, i.e. 

D LL (r) = <{[u(x + r) - u(x)] ■ f } 2 ) = C 2 £ 2 / 3 r 2 / 3 , 
E{k) = C fe e 2 / 3 fc- 5 /3, Ell {k) = C kl e^k~^, ' j 

here -Dll is the second order longitudinal structure function, E{k) is the three dimen- 
sional velocity spectrum, and En(k) is the one dimensional longitudinal velocity spec- 
trum, and k is the wavenumber. The prefactors C2, C k and C k \ in front of the respective 
power laws are the so-called Kolmogorov constants. These constants are widely consid- 
ered universal, i.e. independent of the flow field and the Reynolds number. Once these 
constants are known, they can be used to obtain the energy dissipation rate from the 
measured structure functions or spectra. 

f Email address for correspondence: kxia@phy.cuhk.edu. hk 
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One of the most important parameters in any turbulent flows is the kinetic energy 
dissipation rate. A traditional method for measuring the energy dissipation rate is us- 
ing hot-wire probes while invoking Taylor's frozen flow hypothesis. However, Taylor's 
hypothesis requires that the mean velocity is much larger than the velocity fluctuation. 
Therefore, in turbulent systems with zero-mean velocity or wall-bounded flow with small 
mean velocity, single point velocity measurement could no longer be used to determine 
the dissipation rate. There is a growing consensus recently that the multipoint velocity 
measurement methods, such as particle image velocimetry (PIV) and particle tracking 
velocimetry (PTV), could be used to determine the energy dissipation rate both di- 
rectly and indirectly. The direct method calculates the square of spatial derivatives of 
velocity from the original definition, which is more vulnerable to random errors. The 
indirect method for dissipation rate estimation makes use of the measured structure 
functions or velocity spectra in the inertial range via (|l.lj) . For high Reynolds number 
turbulence the indirect method has great advantages, since the inertial range covers the 
intermediate scales that are easier to measure. Therefore, many experiments with PIV 
or P TV velocity measur ements used this method to determine the energy dissipation 
rate ()Salazar et al. 20081) . As is seen from , in addition to the value of the measured 
Diife), the accuracy of the energy dissipation rate determined from this method also 
depends on the accuracy of Ci u s ed. Th e direct measurements on Ci is very limited. For 
example. lEffinger fc Grossmannl (1987) compared theoretically obtained Ci = 1.7 wit h 
those achieved experimentally Ci = 2.3 ( Anselmet et al. 1984 ; Van Atta &: Chen 1970l ). 
and argued the differ ence is due to overe s timation of dissip ation rate in experiments. In 



many recent studies (jSalazar et all 12008: Vot h et aLl I2002T). th e values of Ci or Ck are 



often obtained from that of Cki using ( Monin fc YaglomT 1975 ) 



C 2 «4.02C fc i, 



55 



(1.2) 



This is because the one-dimensio nal spectra is more accessible experimentally and there 
exists a large set o f data for Cki (ISreenivasanlll995l) . From (|1.2[) and Cki = 0.53 ± 0.055 
( Sreenivasanl 1995 ). one may easily obtain that Ci = 2.13 ± 0.22, which is the most 
widely used value for Ci. A brief survey of the literature shows that (|1.2|) are very often 
used for finite values of R\ when in fact is derived under the idealize d condition of the 
integ ral length L — 00 and the Kolmogorov dissipation scale r\ = ([Monin fc Yaglom 
1975J) , which implies that equation (|1.2[) holds only for an infinite inertial range that may 



be realized only at infinite R\. As shown clearly in (|1.1[) , C2, Ck and Cki are not the 
coefficients of the second-order structure function and spectrum for the whole range of 
scales, but simply the prefactors of their power-law expressions valid only in the inertial 
range. 

In this paper, we examine the relationships among the three Kolmogorov constants 
from the experimental results in different turbulent systems and present two methods 
of determining the i?A-dependent relationships between Ci (Ck) and Cki- The rest of 
this paper is organized as follows. We will first describe the experimental setups in our 
turbulent convection system as well as in other systems in fJ5J The difference between the 
experimental results with previous findings will be discussed in details in H3.ll In H3.21 
we will introduce a method to estimate the limits of finite inertial range and the Ck from 
the second-order structure function. We will extend this method to the estimation of 
the Cki and the ratios among three Kolmogorov constants in £13.31 In addition, another 
method based on the measured inertial range width will be discussed in §3.41 We will 
discuss about our findings and conclude in fJU 
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2. Experimental Setup 

Four sets of data were used in the present study and are labeled as Sets I, II, III and IV 
with their R\ listed in Table I. Set I comes from Rayleigh-B enard turbulent conv ection 
(RBC). Set II comes from an axisymmetric jet experiment ( Anselmet et al. Il984t ). The 
hot-wire for these measurements were positioned at a downstream distance x/d = 25 
from the jet nozzle of diameter d = 12 cm at micro-scale Reynolds number R\ = 536, 
and x/d — 35 at R\ = 852. Set III is from w ind tunnel experiments at NASA Ames 
Research Center ( Saddoughi fc Veeravalli 19941) . The hot wires located on the center line 
of the tunnel ceiling were used to measure the velocity. Set IV is from a direct numerical 
simulation of homogeneou s and isotropic tu rbulence with the resolution ranging from 
128 3 to 1024 3 grid points (|Gotoh et «Lll2002h . 

In the present experiment Lagrangian particle tracking velocimetry (PTV) was used 
to acquire data Set I. It was carried out in a convection cell using water as working fluid. 
The cell is a vertical cylinder with both its height and diameter bein g 19.2 cm. its to p 



and bottom pla tes are made of copper and sidewall made of Plexiglas ISun et al. (2005); 



Ni et al\ ([20121 ). The mean temperature of the bulk fluid is maintained at approximately 
40°C, so the Prandtl number Pr — v/n — 4.4 (v and k being kinematic viscosity and 
thermal diffusivity, respectively). The tracking volume is roughly (5 cm) 3 in the center of 
the cell, which is illuminated by an expanded laser beam. Three fast cameras were used 
to capture the images of seeding particles with diameter 50 /im and density 1.03 g/cm 3 , 
and the acquired images were used as input in a computer program to reconstruct the 3D 
positions of those particles. The Stokes number of the particles in the experiment ranges 
from 10 -3 to 10 -4 , indicating that the particles may be taken as tracers. The resolution 
for the particle positions is ~8 /xm, and the minimum resolvable separation between 
one pair of particles is typically 100 — 200 /im because of the finite particle size and 
diffraction effect. The camera frame rates are 100 fps for all acquired data in Set I, which 
is sufficient to resolve dissipative range properties for the local Kolmogorov time scale 
t v = (i^/e) 1 / 2 , with e being the energy dissipation rate. The velocities were obtained 
by taking direct differentiation of the trajectories, which have been smoothened by a 



are described elsewhere (|Ni et aLll2012l ). 
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Gaussian kernel (IVoth et all 2002). Other details of the setup and the PTV apparatus 



3. Results and discussions 



3.1. Second-order structure function 



Set I was measured in the central region of a cylindrical RBC cell (|Ni et a/.ll2012l) . in 
which the velocity field has been found to be approximately homogeneous and isotropic, 
with the structure functio ns in the inertial range the same as those found for hom ogeneous 
and isotropic turbulence ( Sun et al. 20061 : Zhou et al. l l2008HLohse fc Xiall2010h . Figured] 
shows an example of D^^r) from this Set. The upper inset plots Dll for several values 
of R\, compensated by the dissipative range scaling r 2 , and the lower inset plots Dll 
compensated by the the inertial range scaling r 2 ! 3 . The solid lines in the upper inset give 
the plateau height averaged over data points on the plateau. Using Dr,r, (r) = (e/15v)r 2 
for the dissipative range (r -C 77), we obtain the energy dissipation rate e (|Ni et ctZ.II201 lh . 



It is seen in the lower inset that the plateau that may exist, if at all, must be very short. 
This is due to the limited Reynolds number. Still, we could approximate (jl.l[) to obtain 
C2 by only taking the average over data points near the peak as shown by those solid 
lines. For the data with R\ = 89.7 in figure [TJ this method gives C2 = 1.56, which is 
much lower than the value of 2.13 given by CL2]) with Cfci = 0.53 (|Sreenivasanlll995h . If 
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r / r\ 



Figure 1. A second-order longitudinal structure function (R\ = 89.7). Solid lines: dissipative 
range and inertial range limit. Dashed line: see text. Upper and lower insets: DLL(r) compen- 
sated by r 2 and r 2//3 , respectively. The symbols, from top to bottom, correspond to R\ = 89.7, 
71.8, and 55.0 (Set I). The solid lines are fits to the data in the respective plateaus from which 
e and C'2 are obtained. 



we use Ci — 2.13 to estimate e via (|1.1[) . the obtained value would be 37% smaller than 
that directly measured from dissipative range second-order structure function. This is 
illustrated by the dashed line in figure [TJ which indicates Dll(v) in the dissipative range 
according to e based on C2 = 2.13. 



3.2. Method I 

We now present the first method of determining the i?,\-dependent relationships among 
the Kolmogorov constants from measured structure functions. We note that (|1.2[) were 
derived from the integral transforms between the second-or der structure function a nd 
energy spectra, which, in dimensionless form, are written as (jMonin &: Yaglomlll975l ) 



cos^x sin£x 



m 2 (£r) 3 



(3.1) 



where (3 LL (r/r]) = D LL (r)/v%, ip(r)k) = E(k)/(r)u%), and ipu(r}k) = E u (k) / (r)v%) , with 
Urj = (^e) 1 / 4 the Kolmogorov velocity scale; and x = r/r], £ = krj. One can break up 
the integral in (|3.ip into three parts with respective limits: [0,£i), [^1,^2] and (^2,00), 
where ^1 and ^2 are the yet to be-defined boundaries of the inertial range. The three 
parts correspond to integrations over large-scale, inertial range, and dissipative range, 
respectively. In the idealized case of infinite inertial range, i.e. £1 — > and £2 — ► 00, the 
integrals over the dissipative range and the large-scale vanish and one readily obtains 
(|1.2p by substituting the inertial range results — C22; 2 / 3 , = Cfc£~ 5 ^ 3 and 

< / 3 n(C) = Cfei£ -5 / 3 into (|3.ip and (I3.2p . This proves our assertion earlier that (|1.2p are 
valid only for infinite inertial range, corresponding to infinite R\. 
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Figure 2. (colour online) Circles: compensated Dll for R\ — 89.7 (Set I). Black line: A fit 
of (|3.1[) to the measured structure function. Dashed line and blue line: see text. Inset: More 
examples of compensated Dll with fittings. Blue circles: R\ = 70 (Set IV); green triangles: 
Rx = 600 (Set III). 

To determine £i and £2 for finite R\, we note that the integral over [0,£i) is small 
because the function in the bracket before <p(£) in (|3.1j) approaches zero for very small £ 
(|Monin fc Yaglomlll975l) . The integral over (£ 2 , 00) is also n egligible because the energy 
spectrum in dissipative range decays to zero exponentially (jHeisenber a 1 1 9481) . Thus, to 
a good approximation one may drop these two integrals. The three parameters £1, £2 
and Ck can then be determined by fitting the remaining integral over inertial range 
(with y>(£) = C/c£~ 5 / 3 ) to the measured second-order structure function. The red circles 
in figure [2] shows Dll/ '(re) 2 / 3 (R\ — 89.7, Set I) with the solid line as fitting result, 
which is seen to be excellent. Table I lists the fitted parameters £1, £2 and Ck for all 
data sets. To check whether the fitting result is unique, we changed the fitted value of 
£ 2 to 0.6£ 2 while keeping the other two parameters the same, and used these to evaluate 
the integral in (|3.1[) . which produced the dashed line in the figure. It is seen that the 
right tail of this line is very close to the original fitting line while its left tail shifted 
from the data, suggesting that £2 controls the transition from dissipative scale to inertial 
range. Similarly, we changed £1 to 1.6£i, which gave the blue solid line in the figure. It 
shows that the left tail remains unchanged while the right tail shifted, suggesting that 
£1 controls the transition from the inertial range to the large-scale. It is also obvious 
that Ck controls the plateau value of the compensated second-order structure function 
in the inertial range. Evidently, the measured second-order structure function has three 
different regions: transition zone from dissipative to inertial range, the inertial range, 
and transition zone to large-scale. Thus, for a given compensated second-order structure 
function with left and right tails and a plateau (or a peak), once £1 and £ 2 are fixed, the 
value of Ck is uniquely determined. In the inset of figure [5] we show two more examples, 
which again shows good fitting results. 

3.3. Relations among different Kolmogorov constants 

Although the relationships between C 2 with Ck can be established with a given Dll, it 
is the one-dimensional spectrum that is more experimentally accessible. So it is of more 
practical importance to find relationships between Cki and Ck (C2). For the integral in 
(|3.2p . we take the same treatment for the inertial range and dissipative scales as above, 
i.e. ip(£) — Cfc£ -5 / 3 in [£1, £2] and in (£ 2 , 00). As we are interested in the inertial range, 
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Figure 3. (Colour online) (a) Symbols in the main figure and the inset: ratio of C2/Cki and 
Cfe/Cfci determined from the four data sets with C2 directly from second-order structure func- 
tions, and Ck and Cki from fitting (|3.1|l and (|3.2p to second-order structure functions. Dashed 
and solid lines: see text. The horizontal red line has a value of 4.02. Throughout this figure, 
the symbols are: circles: Set I; triangles: Set II; diamonds: Set III; squares: Set IV. (b) Symbols 
same as in (a) but are plotted as (Cfci/Ca) — 0.25 vs. R\. Dashed line is also the same as in (a). 
Inset: ratio of the parameters a v and ah that characterize the inertial range. Solid line: average 
value of all points excluding diamonds. 



Source 


R\ 




c 2 






6 




Ck 


Cfci 




a v 


I 


55.04 


1 


210 





0515 


0.48 


1 


516 


0.470 


3 


06 


0.48 




71.79 


1 


404 





0332 


0.49 


1 


569 


0.498 


2 


76 


0.49 




89.72 


1 


556 





0278 


0.48 


1 


681 


0.536 


3 


06 


0.48 


II 


536 ± 


2 


067 





0037 


0.68 


1 


804 


0.583 


3 


93 


0.68 




852 


2 


400 





0019 


0.58 


2 


000 


0.647 


3 


51 


0.58 


III 


600 


2 


151 





0013 


0.49 


1 


825 


0.590 


1 


59 


0.49 




1450 


2 


184 





0004 


0.51 


1 


786 


0.578 


1 


50 


0.51 


IV 


70 


1 


893 





0305 


0.38 


2 


180 


0.687 


2 


45 


0.38 




125 


1 


974 





0129 


0.42 


1 


942 


0.625 


2 


16 


0.42 




284 


2 


131 





0045 


0.41 


1 


930 


0.624 


2 


14 


0.41 




380 


2 


138 





0024 


0.40 


1 


856 


0.600 


1 


65 


0.40 




460 


2 


112 





0021 


0.40 


1 


829 


0.591 


1 


83 


0.40 



Table 1. Parameters and fitting result s for the four data set s used in the present work. Set I: 
from this work; Set II: axisymm etric jet (jAnselmet et aZ.II 1984T ) : Set H I: wind tunnel exp eriments 
(Saddoughi & Vccravalli 1994); Set IV: direct numerical simulation ijGotoh et qZ.ll2002l ). 
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kr| 

Figure 4. (Colou r online) Compensated three dimensional spectra taken from Set IV 
<|Gotoh et aZ.ll2002D . From top to bottom, 7? A =460, 380, 284, 125, and 70. The red dashed 
lines are based on parameters from the structure-function-fitting method: they span horizon- 
tally from fi to £2 with height equal to Cfc. The spectra are shifted vertically for clarity. 



we restrict the value of £ in (j3.2[) to £ ^ £1 so we do not need to concern the form of 
<p(£) for £ < £1 [so the obtained spectrum y>n(£) is valid only in the inertial range]. 
With this procedure, we obtain Cki by integrating (|3.2[) with the fitted £1, £2 and Ck, 
and the results are listed in Table I. The values of either C2, or Ck or Cki in the table 
reveal no clear trend with R\ from one set of data to the next or within one set. From 
(jl.ip . we see that the absolute values of Ci and Ck (Cki) determined as the plateaus 
of the compensated structure function/spectra will be affected by the errors in e, but 
such errors will be canceled out if we take their ratios. As shown by the symbols in 
figure |^a) and the inset, this appears to be the case, i.e. clear trends emerge with R\ 
for these ratios. To find an empirical relation between C2 and Cki, we fit the symbols in 
figure |3][a) with C^/Cfci = 1/(0.25 + AR^J. The result is shown as the dashed-line in the 
figure, with the fitting parameters A = 1.95 and j3 = —0.68. In figure [3jb) we plot the 
data as (Cki/Cz) — 0.25 vs. R\ with the dashed-line representing the fitting, which shows 
that the plotted quantity indeed follows a power-law within the data range, suggesting 
the fitting curve in figure EJa) can describe well the behavior of CifCki- 



3.4. Method II 

We now introduce a second method, based on experimental values of L/rj, to determine 
the ratios of CijCki and Ck/Cki, but not the individual constants. Let £1 = aL(r)/L) 
and £2 — a n , respectively. Table I lists the values of and a n for the four sets of data. It 
is clear that these numbers depend on R\, but their ratio appears to be independent of 
R\ [see the inset of figure [S^b)] (the diamonds are from Set III a nd their large deviations 



from the others may be due to the noise in the original data (jSaddoughi h Veeravalli 
19941) . The horizontal line in the figure represents a v /aL — 0.183 and is an average based 



on Sets I, II and IV. [If all data points are included, the value will be 0.21 instead. 
But the difference in the obtained CijCki in using either one of these values is small.] 
It thus appears that although both the inertial range itself (characterized by L/rj) and 
its start and end points (represented by £1 and £2 in fc-space) depend on R\, the ratio 
a n /aL — (£2 / £1) / (L / rf) seems to be universal. We note that if £1 and £2 are known for 
a given R\, we can then simply substitute (jl.ip into the integrals in (|3.1[) and (I3.2j) to 
obtained the ratios, without the need for either structure function or spectra data. For 
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those data collected by Cleve et al. ( 20041 ). we obtain values of L/rj, which are based on 



a number of experiments with R\ varying from 85 to 17,000. A power-law fit to these 
gives L/rj = 0.37 R\ 27 . It is seen from Table 1 that a v is generally between 0.4 and 0.6. 
Taking a v — 0.5, and with a v /aL = 0.183, we obtain £2 and £1 as functions of R\. For 
each of these pairs, we obtain C^/Cfc and Ck/Cki via equations (jl.ip . (|3.1[) and (13. 2[) . 
and therefore C^jCki- The results are shown as the solid lines in figure [31(a) and the 
inset, respectively. One sees there is excellent agreement with those obtained from fitting 
structure function data, as well as with the fitted curve. We remark that if one takes a v to 
be any value between 0.4 and 0.6, the obtained CijCki is essentially the same, especially 
when R\ > 100. It should be note that in most real flows there exist no sharp boundaries 
between the three regions, i.e. the dissipative range, the inertial range and the large 
scales, but only transition zones. In this respect, the parameters £1 and £2 introduced 
here should be viewed as a mean to obtain the ratio of the Kolmogorov constants, which 
we have shown to be largely determined by the ratio £2/^1 and is insensitive to the 
specific values of these parameters if they are within a range. 

As a self-consistent check of the structure-function-fitting method, we compare the 
obtained £ 1; £ 2 and Ck with 3D spectra, which are from the same DNS study as the 



structure functions from Set IV but were obtained separately ([Gotoh et ali 120021 ). Fig- 
ure [J] plots 5 compensated spectra with R\ from 70 — 460. The horizontal dashed lines 
are determined by the three fitted parameters, i.e. they start at £1, end at £2, and have 
heights equal to Ck ■ It is seen that the method can predict accurately not only the height 
of the spectrum plateau, but also the position and the width of the inertial range except 
the ending points, which are a bit larger than those shown in the spectra. 



4. Discussions and Conclusion 

In the paper we have introduced two different methods that give the relationships 
among the three Kolmogorov constants C2, Ck, and Cki- The first method can be used 
to determine Ck and Cki from the experimentally measured second-order structure func- 
tions by splitting the integral relation between the structure function and the energy spec- 
trum into three different regimes: dissipative range, inertial range, and range for large 
scales. It is found that the ratio C^jCki exhibits clear dependence on the micro-scale 
Reynolds number R x . The ratio could be well fitted as (C k i/C 2 - 0.25) = 1.95i?^ 68 . 
The second method directly determines the ratios among the three Kolmogorov con- 
stants from the experimentally-determined width of the inertial range without fitting 
to the second-order structure function. Our results reveal that the widely-used relation 
C2 = 4.02Cfci holds only asymptotically when R\ > 10 5 . 

It is generally known that C2, C'k and Cki are asymptotic values, which will become 
a constant for high enough Reynolds number. Since the Cki is generally considered as 
constant for R\ > 50 and it is more experimentally available, C2 are usually obtained 
from C2 — 4.02Cfci for R\ from hundreds to thousands. However, this is only valid if C2 
equals to 4.02Cfei for all R\. And of course, the C2 should be constant for R\ > 50. In 
this work, we have shown that the C^jCki are Reynolds number dependence at R\ as 
high as 10 4 , which indicates that C2 does not reach asymptote at least for R\ < 10 4 . 

The reason for the stronger R\ dependency of C2 is because there is no clear inertial 
range for second-order structure function for R\ even at 10 4 . It is found that if the 
structure function showed in log-log plot, there may have clear inertial range, however 
if the SF2 is compensated and plotted in log-linear scale, the inertial range may not 
be that apparent. It is consistent with previous findings that there is no inertial range 
for SF2 at R\ — 19500 with a more sensitive local slope test. In all R\ examined, we 
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approximated to obtain C'2 by only taking the average over data points near the peak 
in the inertial-range-compensated structure function. The reason we could still use this is 
because, physically, the inertial range does exist for all R\ we used. There is an inertial 
range in three-dimensional energy spectrum which is the prerequisite for our method. 
But after transformation through (13.11) . the range is compressed in structure function 
form. And there are two benefits for this method. First, it could help us to determine 
when the C2 will reach the asymptote. By extrapolation from our works, we estimate it is 
roughly at R\ = 10 5 . Second, even with this approximated C2, the experimentalist still 
could obtain the dissipation rate through the structure functions. The relations produced 
by our analysis would give a more accurate result than simply using C2 = 4.02Cfci for 
estimation. 

It is foun d before that both C 2 and Ck i have great uncert ainties such as those com- 
piled by M onin fc Yagloml (|l975l ) and by ISreenivasanl (|l995f) . The uncertainties would 
come from the method in determining dissipation rate or the energy injection method in 
different turbulent systems. We may ascribe both ways to the error of energy dissipation 
rate. When comparing the Kolmogorov constants obtained in different systems, if there 
is a large uncertainty, any trend below would be hidden in the scatter exhibited by the 
data. Therefore, in this work, we consider the Reynolds number dependencies for the 
ratios rather than the value for each constant separately, since the contributions from 
energy dissipation error will be cancelled with each other. 

We gratefully acknowledge support of this work by the Research Grants Council of 
Hong Kong SAR (No. CUHK404409 and N_CUHK462/11). 
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